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ABSTRACT 

The weak lensing signal (cosmic shear) has been shown to be strongly contaminated by 
the various types of galaxy intrinsic alignment (IA) correlations, which poses a barrier 
to precision weak lensing measurements. The redshift dependence of the I A signal 
has been used at the 2-point level to reduce this contamination by only measuring 
cross-correlations between large redshift bins, which significantly reduces the galaxy 
intrinsic ellipticity - intrinsic ellipticity (II) correlation. A self-calibration technique 
based on the redshift dependences of the IA correlations has also been proposed as 
a means to remove the 2-point IA contamination from the lensing signal. We explore 
here the redshift dependences of the IA and lensing bispectra in order to propose a 
self-calibration of the IA auto-correlations at the 3-point level (i.e. GGI, Gil, and III), 
which can be well understood without the assumption of any particular IA model. 
We find that future weak lensing surveys will be able to measure the distinctive IA 
redshift dependence over ranges of |Az p | ^ 0.2. Using conservative estimates of photo- 
z accuracy, we describe the 3-point self-calibration technique for the total IA signal, 
which can be accomplished through lensing tomography of photo-z bin size ~ 0.01. 
We find that the 3-point self-calibration can function at the accuracy of the 2-point 
technique with modest constraints in redshift separation. This allows the 3-point IA 
auto-correlation self-calibration technique proposed here to significantly reduce the 
contamination of the IA contamination to the weak lensing bispectrum. 

Key words: gravitational lensing - cosmology 



. !_H . 1 INTRODUCTION 

r j_] ■ Originally just a prediction of general relativity, gravitational lensing has recently emerged as an independent means to make 
C^ , precision astrophysical and cosmological measurements, which are blind to the exact nature of the lens mass. This is particu- 
larly true for weak lensing due to large scale structure (cosmic shear), which is able to map structure composed of both visible 
matter as well as dark matter, which is primarily detectable through its gravitational signal only. The importance of this new 
probe has spurred the development of a new generation of ground- and space-based surveys suited for precision weak lensing 



measurements. These ongoing, future, and proposed surveys (e.g. CFHTLS (http://www.cfht.hawaii.edu/Science/CFHLS/), 
DES (http://www.darkenergysurvey.org/), EUCLID (http://sci.esa.int/euclid/) , HSC (http://www.naoj.org/Projects/HSC/), 
HST (http://www.stsci.edu/hst/), JWST (http://www.jwst.nasa.gov/), LSST (http://www.lsst.org/lsst/) , Pan-STARRS 



(http://pan-starrs.ifa.hawaii.edu/), and WFIRST (http://wfirst.gsfc.nasa.gov/)) promise to provide greatly improved mea- 
surements of cosmic shear using the shapes of up to billions of galaxies. These cosmic measurements allow us to not only 
characterize the equation of state of dark energy, but when combined with other probes can improve constraints on the 
equation of state of dark ene r gy an d the matter fluctuation amplitude parameter by fa c tors o f 2 to 4 (see for example 



Zaldarriaga, Spergel fc Seliak I hmij): lEisenstein. Hu fc Tegmark I Jl999h; feu fc Tegmark I Jl999lh iBacon Refregier fc Ellis 



f ldarriaga, Spergel fc Seliak I Jl99 , i]); 
IQOlk lVan Waerbeke et al. I (^OOol hTR 



2000t l: IVan Waerbeke et al. l(|200(]h : lRhodes. Refregier fc Groth I (|200lh : lHoekstra et al. l|2002h :lHu l |2002J h lVan Waerbeke et al 
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Figure 1. The 2-point intrinsic alignment correlations. Galaxies which are intrinsically aligned are coloured in blue and labeled I, while 
galaxies which are lensed are coloured in red and labeled G. The lower right panels represent the view of the system on the sky, while 
each panel preceding it is at some distinct redshift where Zj < Zj . If the two galaxies are spatially close, at nearly the same redshift 
and angular position on the sky, they can be aligned by the tidal force field of the same nearby matter structure (labeled DM in the 
figure). This is shown as the II correlation. If instead a matter structure causes both the alignment of a nearby galaxy and contributes 
to the lensing signal of a background galaxy, this produces an anti-correlation with negative sign between the cosmic shear and intrinsic 
ellipticity, since the tidal force and gravitational lensing tend to align the galaxy shapes in orthogonal directions, and is shown as the GI 
correlation. 



2002) 



2005) 



Brown et al~l d2003ri:ljarvis et al. I d2003h:IPen et al. N2003I): llshak N2005l);lMassev et al N2005I ) ; lUpadhve. Ishak fc Steinhardt 



Ishak 



( 20071 ): iFu. Wu fc Yul(|2009l ) lJoudaki. Coorav fc Holz I J2OO9F 



Schrabback et al. 



1 2OIO1 ) and references therein.) 



Weak lensing has also been shown t o be very useful to test the nature of gravity at cosmological distance s c ales (see for 
examp l e the partial list Song ( 2005 ); Capozziello. Cardone fc Troisi H2006I) ; llshak. Upadhve fc Spergel ( 2006 ); Zhao et al 



j2006f);lHuterer fc Linder I d2007l);lLinder fc Cahn I d2007TI;IZhang et al. I dgOOThjAcq 






I'M 

18 s ! 



uaviva et al. I (120081) ;lDaniel et al. I (120081); 
: I J2009I); IZhao et al. I J2009I); bean fc Tangmatitham ' 



( 2Q10l) ; lDaniel et al. I J2010l ); lTereno. Semboloni fc Schrabback I (J2010T l; lDossett. Moldenhauer fc Ishak I (J201ll ).) 

The 3-point cosmic shear correlation and the shear bispectrum have been shown to break degeneracies in the cosmological 
parameters in addition to the constraints obtained from the 2-point cosmic shear correlation and the corresponding shear 
p ower spectrum that the power spectrum alone does not ( Takada fc Jain 1120031 ; IVafaei et al. II2OIO1 ). For example, the results 
of iTakada fc Jain I (|2004l ) showed that a deep lensing survey should be able to improve the constraints on th e dark energy pa- 



ramet ers and the matter fluctuation amplitude by a further factor of 2-3 using the bispectrum. Most recently. ISemboloni et al 



1 20101 ) derived parameter constraints by measuring the third order moment of the aperture mass measure using weak lensing 
data from the HST COSMOS survey. They found independent results consistent with WMAP7 best-fit cosmology and an 
improved constraint when combined with the 2-point correlation. In addition to improved parameter constraints, by defini- 
tion the bispectrum also allows us to explore information about non-Gaussianity in t he universe that is inaccessible at th e 
2-point level , prov i ding constraints on the degree of non-Gaussianity (see for example M^jan^s^ i _yenJB_fc_JjrnBnez_ (2000) 



Verde et al.l(|200ll ); lTakada fc Jainl J2004h ; [jeong fc Komatsu1(|2009l ); lrluterer, Komatsu fc Shandera I (|2oiol ) 
( 201 ll ) and references therein.) 



Munshi et al. 



As we describe below, we extend in this paper the 2-point self calibration technique proposed by IZhang I (|2010bl ) to the 
3-point intrinsic alignment auto-correlation bispe ctra between ga laxie s in a single red s hift bi n. This technique is different 
from the cross-correlation techniques proposed in IZhang I (|2010al ) and I Troxel fc Ishak I (|2012al ) , instead using differences in 
the redshift dependencies of the intrinsic alignment and lensing bispectra to self-calibrate the intrinsic alignment signal. These 
2-and 3-point intrinsic alignment correlations constitute a contaminant to the lensing signal and must be isolated and removed 
to avoid biasing the cosmological information contained within the cosmic shear power spectra and bispectra. 

Cosmic shear measurements are in fact limited in precision by several systematic effects which must be ac c ounte d 
for in order to make full use of the potential of future weak lensing surveys (see for example ICroft fc Metzler I (12 000) 



Heavens. Refregier fc Hevmans N2000l);lBacon et al N200ll);ICatelan. Kamionkowski fc Blandford ll200ll) ; lErben et al N200ll ) 



Van Waerbeke fc Mellier I (|2003l ); lHevmans et al. I (|2004D ; llshak et al. l(|2004l ); lTakada fc White I (|2004l ) and references therein) 



Bernstein fc Jarvis I d2002Tl; iBrown et al. I J2OO2I) ; iKing fc Schneider" d2002rTlIirata fc Seliakl d2003al ); iRefregier I "(J2003I ) 
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Figure 2. The 3-point intrinsic alignment correlations. Galaxies which are intrinsically aligned are colored in blue and labeled I, while 
galaxies which are lensed are colored in red and labeled G. The lower right panels represent the view of the system on the sky, while 
each panel preceding it is at some distinct redshift where zt < Zj < z^ . The III correlation is between the intrinsic ellipticities of three 
spatially close galaxies which arc intrinsically aligned by a nearby matter structure (labeled DM in the figure). If instead two spatially 
close galaxies are intrinsically aligned by a nearby matter structure which contributes to the lensing of a third galaxy in the background, 
we label this the Gil correlation. Finally, the GGI correlation is where two galaxies are lensed by a structure which intrinsically aligns 
a third galaxy in the foreground. Unlike the 2-point correlations, the sign of the GGI and Gil correlations can depend both on triangle 
shape and scale. 

and one of the serious systematic ef fects of weak lensing is t h is correlated intrinsic alignment of galaxy ell i pticities, which act 



as a n u isance factor (see for e xampl e Croft fc Metzler (120001); Catelan. Kamionkowski fc Blandford 



)7) 



j200ll);lBrown et al I JboOgh ; I Jing I hOQJ) ; king fc Schneider I (|2002ri;lHevmans fc Heavens I (hOQot); 



Ci 

<; 

)5) 



King fc Schneider I J2003h:|Hirata fc Seliak N2004ri;lKing N2005h:lrlevmans et al.' J200d):lMandelbaum et al. N200rj);lBridle fc King 



(2001); C rittenden et al . 
Hirata fc Seliak d2003bF: 



Hirata et al. (120071); Semboloni et al. (hoosf); Faltenbacher et al. ( 2006 ); Okumura fc Jing I J2009h ; Joachimi fc Bridle 



2007) 



2010) 



Joachimi et al. 



l|2010l) ; lKirk. Bridle fc Schneider I J2010l);lBlazek. McQuinn fc Seliak I (hoilhjKrause fc Hirata I J201J) ; 
Troxel fc Ishak I ( 2012bl ) and references therein). For example. iBridle fc King I ( 20071 ) and I Joachimi fc Bridle I ( 20101 ) showed 
that if intrins i c alig nment is ignored, the determination of the dark energy equation of state is biased by as much as 50%. 



Hirata et al. 



(|2007f ) found that the matter power spectrum amplitude can be affected by intrinsic alignment by up to 30%, 
showing the importance of developing methods to isolate the intrinsic alignment and remove it from the cosmic shear signal. 
There are two 2-point intrinsic alignment correlations. The first is a correlation between the intrinsic ellipticity of two 
galaxies, known as the II correlation. If the two galaxies are spatially close, they can be aligned by the tidal force field of 
the same nearby matte r structure. The second intrinsic alignment correlation, known as the GI correlation, was identified by 
Hirata fc Seliak I (|2004f ) and is due to a matter structure both causing the alignment of a nearby galaxy and contributing to the 
lensing signal of a background galaxy. This produces an anti-correlation between the cosmic shear and intrinsic ellipticity, since 
the tidal force and gravitational lensing tend to align the galaxy shapes in orthogonal directions. We represent these correlations 
diagrammatically in Fig. [T] The GI correlation has been measured in various subsets of the SDSS spectr oscopic and imaging 
sampl es by vari ous groups. A detect ion of the large-scale GI correlation in the SDSS was reported by iMandelbaum et al. I 
(2006!), and then lHirata et al. I ( 20071 ) found an even stronger GI correlation for Luminous Red Galaxies (LRGs). It was shown 
in these papers that this contamination can affect the lensing measurement and cosmology up to the 10% level and up to 30% 
in some cases for the matter fluc tuation amplitude. Th is finding was confirmed by numerical simulations, where a level of 



contamination of 10% was found ( Hevmans 



et al. 



200rj). 



Further measurements of the GI correlation were made in the SDSS 



dataset by iFaltenbacher et al. I ( 20091 ) and lOkumura fc Jing I ( 20091 ) . Most recently, Ijoachimi et al. I ( 20101 ) measured strong 



2-point intrinsic alignment correlations in various SDSS and MegaZ-LRG samples. 

In a similar way, when we consider three galaxies and the related 3-point correlation, the cosmic shear signal (GGG 
bispectrum) also suffers from contamination by the 3-point intrinsic alignment correlations. The first is the III correlation 
between intrinsic ellipticities of three spatially close galaxies which are intrinsically aligned by a nearby matter structure. The 
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second is the Gil correlation, where two spatially close galaxies are intrinsically aligned by a nearby matter structure which 
contributes to the lensing of a third galaxy in the background. Finally, there is the GGI correlation, where two galaxies are 
lensed by a structure which intrinsically aligns a third galaxy in the foreground. Unlike the 2-point correlations, the sign of the 
GGI and Gil correlations d epend both on triangle s hape and scale. The 3-point intrinsic alignment correlations are represented 
diagrammatically in Fig. l2l lSemboloni et al. I (2003) showed that lensing bispectrum measurements are typically more strongly 
contaminated by intrinsic alignment compared to the lensing spectrum measurements, and that the contamination from the 
3-point intrinsic alignment correlation can be as large as 15 — 20% compared to the GGG lensing signal. The 3-point intrinsic 
alignment measurements are not only useful for constraining their contamination to 3-point lensing measurements, but are also 
useful for constraining models of intrinsic alignments and therefore constraining the contamination to all lensing measurements 
(including 2-point correlations) which will dominate the science cases of upcoming surveys. 

While the II, III, and Gil intrinsic alignment correlatio ns can be g reatly reduced with photo-z's by using cross-spectra 
of galaxies in two different redshift bins (see for example iRefregier I ( 20031 )) so that the galaxies are separated by large 
enough distances to assure that the tidal effect is weak, this does not work for the GI and GGI correlations which remain 
strong between galaxies at different redshifts and large separations. The GI correlation and methods for its removal have 
been the to pic of several recent sci e ntific publica tions and we review these briefly. Initially, s ome first suggestions were 



discussed bv lHirata fe Seliak I ( 2004 ). iKing I (2005 ) extended the a ppro ach of template fitt i ng by iKing fe Schneider I (2003) 



to include a treatment of the GI correlation. iBridle fc King I ( 20071 ) and I Joachimi fc Bridle I ( 20101 ) investigated the effects of 



the GI correlation on cosmological parameter constraints by assuming a model of the GI intrinsic alignment that is binned 



in redshift and angular frequency with some free parameters that are marginalized over. iKirk. Bridle fc Schneider I (2010J) 



performed a cosmological constraint analysis where modeling of intrins ic alignment was included, sh o wing a significant effect 
on the amp litude of matter fluctuat i ons. U sing a geometrical approach. I Joachimi fc Schneider 1 1 20081 ) . I Joachimi fc Schneider 



1 2009 ). and I Joachimi fc Schneider I ( 20101 ) proposed a nulling technique to remove the GI intrinsic alignment contribution by 
exploiting the redshift dependence of the correlations, but it was found that the technique throws out some of the valuable 
lensing signal. Most recently, the nulling tech nique has been applied at the 3-po int level f or the GGI corr elation, but again with 
similar signal loss to that at the 2-point level ( Shi. Joachimi fc Schneider 1120101 ). Finally. IZhang I (2010aj) proposed a technique 
to self-calibrate the GI intrinsic alignment signal by using additional galaxy density (cross-)correlations w hich are already 



present in weak lensing survey measurements. Th is approach was suc c essful ly extended in recent work by iTroxel fc Ishak 



1 2012al ) to the GGI cross-correlation bispectrum. Ijoachimi fc Bridle I (2010j) applied an approach like the self-calibration, 
using correlations between lensing, intrinsic alignment, number density, and magnification effects to constrain cosmological 
parameters. They found that the extra information from the additional correlations can make up for the additional free 
parameters in the i ntrinsic alignme nt so that the contamination can be removed without loss of constraining power. 

Most recently, IZhang I (J2010bJ) showed that redshift dependencies of the intr i nsic ali gnment spectra can allow further 
improvements to the calculation of the intrinsic alignment contamination. IZhang I (2010bJ) demonstrates the strong redshift 
separation dependence of the 2-point GI and II intrinsic alignment signals by considering the lensing and intrinsic alignment 
spectra due to galaxies in a single redshift bin that are at fixed photo-z separations Az p relative to some mean redshift. This 
results in a relative change in magnitude of the GI and II spectra of 50-60%, but only a few percent for GG at separation 
Az p = 0.2. This corresponds to a decrease in 10% of the total ellipticity power spectrum, which is identifiable in proposed 
surveys. Parameterizing the intrinsic alignment spectra as proportional to various galaxy density (cross-)spectra, this allows 
us to calculate and remove the intrinsic alignment component of the lensing signal. We propose in this work a method to 
extend this self-calibration technique to the 3-point statistics, using instead the redshift dependencies of the GGI, Gil, and 
III bispectra, in order to calculate and remove the intrinsic alignment contamination to the lensing bispectrum. 

We organize the paper as follows. In Sec. [2] we briefly discuss the necessary survey parameters and lensing formalism 
for the bispectrum. In Sec. 12.11 we discuss several methods to vary the redshift separation of galaxies in the bispectrum 
and compare the resulting redshift dependencies. Section T2.2I describes the analytical motivation and framework for the self- 
calibration using the redshift dependencies of the GGG, GGI, Gil, and III bispectra. We explore the performance of the 
self-calibration in Sec. 13.11 and discuss possible additional sources of error in the analysis in Sec. 13.21 Finally, we summarize 
the self-calibration and its applicability to future weak lensing surveys in Sec. [3] and discuss its impact in Sec. [4] 



2 REDSHIFT SEPARATION DEPENDENCE IN THE BISPECTRUM 

The galaxy intrinsic alignment and cosmic shear bispectra have very distinct dependences on the redshift separation and 
orientation of the three galaxies, which we will explore and quantify in order to identify and measure the contribution of the 
intrinsic alignment to the measured ellipticity bispectrum. We work with the bispectrum B a ^ J (£; z p , z p , z p ), where a,/3,j 
6 G,I,g measures a at photometric redshift z p , /? at z P , and 7 at z P . We have denoted the lensing convergence as G, the 
galaxy intrinsic ellipticity as /, and the galaxy density as g. Since we are only interested in exploring the redshift separation 
and orientation dependence of the bispectra, we will work at a fixed multipole I and mean redshift z p — (z p + z p + z p )/i. 
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In order to best compare to work on the power spectrum by I Zhang I (|2010br ). we will use equilateral triangles with £ = 1000 
and z p — 1.0, unless otherwise stated, for which values the bispectrum will be measured to high confidence in planned weak 
lensing surveys. Unlike for the power spectrum, there are several possible choices in exploring the dependence on redshift 
separation, which we will discuss and compare in section 12.11 

We assume in our calculations a standard, flat ACDM universe. In the Born approximation, the convergence k of a source 
galaxy at comoving distance X G and direction 9 is then related to the matter density S through the lensing kernel Wl(z' , z) 
by 

K0) = f XG S( XL , 6)W L {xL,XG)d X L. (1) 



The 3D matter bispectrum is then defined from the convergence as 



{k(£i)k(£ 2 )k(£ 3 )) = (2n) 2 8 D (£ 1 + £ 2 +£ a )B(e 1 ,e 2 ,e 3 ), 
where (• ■ • } denotes the ensemble average and S D (£) is the Dirac delta function. For the bispectrum, 8 L 



(2) 
■f £2 + £3) enforces 

the condition that the three vectors form a triangle in Fourier space. Under the Limber approximation, we can express the 
2D angular auto-correlation bispectrum as 

rx W a ^( x ';xi,X2,X3) 



B' 



*P~t 



(£; Zi , z 2 , 23 ) — 



-B al3l {k;x)dx , 



(3) 



/o X 

where for example when a = fi — 7 = G, BcGG(k;x') i s the 3D matter bispectrum shown in Eq. [2] However, generally 
a, P, 7 £ G, I, g, where the additional intrinsic alignment (I) and galaxy (g) bispectra are calculated as described further in 
Sec. 12.11 The redshift is simply related to \ through the Hubble parameter, H(z). We can then write the weighting function 
in terms of redshift as 



W a ^(z( X );z[( Xl ),z[( X2 ),z[( X3 )) 



W a (z,z P )W^(z,z P )W^(z,z P ), 



W G (z,z p ) 



W L (z',z)p(z'\z p )dz', 



(4) 

(5) 

W 9 (z,z P )=p(z\z P ), (6) 

where p(z\z p ) is the photo-z probability redshift distribution (PDF). In order to quant itatively examine t he red shift separation 
dependence of the bispectra, we assume a specific form of the PDF modeled after iMa fc Bernstein I ( 20081 ) with which to 
describe the photo-z uncertainty in a future weak lensing survey 



W'{z 



* P ) 



P{z\z 



1 — Pc 



2-k<j(z p ) 



exp 



(*■ 



2a 2 (V 



+ 



Peat 



2iva(z p ) 



exp 



(z — fbizsz p y 



2a 2 



(7) 



The fraction of outlier galaxies is given by p ca ,t , with a true redshift which is biased by a factor /bias- We adopt similar values to 
Zhang, with <j(z p ) = 0.05(1 + 2 F ), p ca t = 0.02, and / b i as = 0.5. These parameters are chosen to represent the expected photo-z 
accuracy of a stage IV weak lensing survey at z p ~ 1. We maintain the very conservative fraction of catastrophic outliers 
used by Zhang, which is a factor of 20 greater th an the required fraction in a stage IV dark energy survey, where p C at < 0.1% 
( Bernstein fc Huterer 1 12009 ; iHearin et al. luOlOT ). Typically we assume that B (£; 21,22,23) is zero for 21,22 < 23 and 
that B (£; 21, 22, 23) is zero for 21 < 22, 23 or 22 7^ 23, due to lensing g eometry and the redshi ft separation and orientation 
dependence of the intrinsic alignment signal discussed in detail in Sec. l2.1l ( Troxel fc Ishak Il2012al ). However, due to sometimes 
large photo-z error, this is not always the case, so when referencing the GGI and Gil correlations, we will instead work with 
the sums B GGI + B GIG + B IGG and B Gn + B IGI + B IIG , respectively. 



2.1 Evaluating the redshift separation dependence 

We calculate th e requir ed le nsing, intrinsic alignmen t, and galaxy bispectra using the relations and methods developed in 
Troxel fc Ishak I (|2012al ) and lTroxel fc Ishak I (J2012bJ). We ass ume a deterministic galaxy bias for the galaxy bispectra, while 
the intrinsic alignment signal is calculated using the model of ISchneider fc Bridle I ( 2009T ) (SB09), which is based on the halo 
model prescription. We have used th e fiducial parameters of their fitting formulae as listed in their Tables I & II, with G\ 
es timated by comparison to Fig. 2 of iHirata fc Seliak I (J200J). By design, this model reduces to the linear alignment model 
of iHirata fc Seliak I (J200J) at large scale, but aim s for a more m otivated model of intrinsic alignment at small scales. For 
comparison, we also use the toy model adopted bv lZhang I (J2010bJ), where the intrinsic alignment spectrum has a simple bias 
b 1 (k, 2) oc (1 + A 2 n (fc, z)) a (a £ [0, 1/2]), where A 2 n (k,z) is the three-dimensional matter power spectrum. This allows us to 
evaluate the behavior of the redshift separation and orientation dependence of the bispectrum for a wide range of intrinsic 
alignment dependence on both scale and redshift. 

We relate the flat-sky bispectrum to the all-sky bispectrum through the Wigner-3j symbol, where 



B 



a/37 



(2 





(2ll + l)(2l2 + l)(2l 3 + l) pa) 3 7 



4tt 



B^(£ 1 ,£ 2 ,£ 3 ). 



(8) 
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Method 1 


Method 2 


Method 3 


Az p 


0.1 


0.2 


0.1 


0.2 


0.1 0.2 


GGG 


0.5% 


1.5% 


0.5% 


1.5% 


0.5% 2% 


GGI 


12% 


48% 


12% 


37% 


20% 89% 


Gil 


1% 


6% 


5% 


32% 


8% 42% 


111 


23% 


64% 


28% 


72% 


30% 77% 



Table 1. The Az p dependence of B GGG (GGG), B GGI + B GIG + B IGG (GGI), B Gn + B IGI + B IIG (Gil), and B UI (III) for the 
three methods of varying rcdshift separation described in Sec. 12. ll The values listed are the ratio \B(Az p )/B(Az p = 0)|, showing the 
relative change in magnitude between Az p = 0.0 and Az p = 0.1 and between Az p = 0.0 and Az p = 0.2. Method 1 shows a similar 
Az p dependency to the power spectrum, while Method 3 clearly displays the most distinct Az p dependency. 



We calculate the Wigner-3j symbol following the approximation given in Appendix A of iTakada fc Jain I ([2004]). We then 



compute t he three-dimensional bispectrum due to nonlinear gravitational clustering, Bsjki, k 2 , k 3 \ y), following the fittin g 



formula of lScoccimarro fc Couchman I ( 2001m with coefficients F 2 (fci, k 2 ) described in Section 2.4.3 of ITakada fc Jain I (2004) 



B s (ki,k2,k 3 ;x) = 2Ff(k ly k 2 )Ps(k 1 ; X )Ps(k2;x) + 2 perm. (9) 

In order to approximate the three-dimensional intrinsic alignment bispectra, we make a direct expansion of this method, using 
the intrinsic alignment power spectra instead of the nonlinear matter power spectrum, where 

B S5¥ (ki,k2,k 3 ;x) = 2Ff{k 1 ,k 2 )P s¥ {k 1 -x)P6{k 2 -x)+2Ff{k 2 ,h i )P & {k 2 - x )P 5¥ {k 3 -x) (10) 

+2Ff(k 3 , kx)P s¥ (k 3 ; X)P S¥ (foiX) 
B s¥sjI (k 1 ,k 2 ,k 3 ;x) = 2Ff(k 1 ,k 2 )P ¥ (k 1 ;x)P S¥ (k 2 ;x) + 2Ff(k 2 ,k 3 )P s¥ (k 2 ;x)Ps ¥ (k 3 ;x) (11) 

+2Ff(k 3 ,k 1 )P ¥ (Asa; x )Ps-- l i (ki;x) 
B ¥ (ki,k2,k 3 ;x) = 2Ff(k 1 ,k 2 )P ¥ (k 1 - X )P ¥ {k 2 ;x) + 2 perm. (12) 

We find that this treatment gives reasonable results for the intrinsic alignment bispectra. 

We explore several methods of varying the redshift separation and orientation in the bispectrum. We can define the 
bispectrum triangle configuration by the redshift of its vertices (z p , z 2 ', and z P ) and the £- modes corresponding to the 
angular separation of its sides. From these values, we can derive the mean redshift (z p ) and the side lengths measured in 
redshift of each triangle (Az p , Az 2 , and Az P ), where the side Az p is defined opposite the vertex n with photo-z z p . In all 
cases, we keep z p constant. The first (Method 1) is most similar to the method employed by Zhang for the power spectrum. 
We set z p = z p constant, and instead vary z p and z P such that Az p = Az P = Az 2 /2. The quantity Az p is then the 
separation that we vary, which completely defines the redshift configuration of the triangle in this method. 

In Methods 2 and 3, we set z p = z 2 such that Az P — 0. We then vary Az p — Az p , which completely defines the 
triangle configuration. In the first case, Method 2, we vary the separation Az p = Az 2 with z p increasing in redshift. In the 
second case, Method 3, we instead vary the separation Az p = Az p with z P decreasing in redshift. We can also explore more 
complicated methods where Az P is constant and non-zero. We then vary Az 2 = \Az P — Az p |. These more complicated cases 
do still show a large change in the intrinsic alignment bispectra with varying Az 2 , but are less useful for a self-calibration of 
the intrinsic alignment signal as there is not a simple relationship between the intrinsic alignment and galaxy bispectra as we 
discuss in Sec. [3] 

We summarize the Az p dependence of Methods 1-3 in Table 1, showing the percent change in the magnitude of the 
bispectra from Az p - 0.0 to Az p =0.1 and from Az p = 0.0 to Az p = 0.2. With Method 1 we find a resulting Az p 
dependence in the GGI and III bispectra which is similar to that found by Zhang for the GI and II power spectra, though 
with almost no detectable dependence in the Gil bispectrum. This similarity and the reduced dependence in Method 1 when 
compared to Methods 2 and 3, which are described below, can be understood from the way we have varied the redshifts in 
this case, keeping one vertex at a fixed redshift z p , which is very similar to the 2-point case where only two points are varying 
in redshift. Both Methods 2 and 3 show a noticeable increase in Az p dependence, though the GGI bispectrum in Method 
2 is actually less dependent than in Method 1. All three demonstrate a nearly identical percent change over a separation of 
Az p = 0.2 in the GGG bispectrum, at the 1-2% level. 

Method 3 clearly has the strongest Az p dependence. We find this is significantly stronger than for the 2-point case and 
both the previous Methods, and so we will focus our discussion on Method 3 as the best candidate for a self-calibration of 
the intrinsic alignment signal. As described in Table 1, we find a relative change of 8, 20, and 30% in the Gil, GGI, and III 
bispectra, respectively, between Az p — 0.0 and Az p = 0.1. However, this increases greatly such that between Az p = 0.0 and 
Az p = 0.2, we find a very large relative change of 42, 77, and 89% in the Gil, III, and GGI bispectra, respectively. The change 
in the GGG bispectrum is only 2% between Az p = 0.0 and Az p = 0.2, by comparison. This very strong Az p dependence will 
allow us to develop a means to self-calibration of the intrinsic alignment signal, which we describe in the following sections. 
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Figure 3. The weighting functions W a ^ evaluated at I = 1000 and z p = 1.0. The lines labeled GGG correspond to W GGG , GGI 
to the sum W GGI + W GIG + W IGG , Gil to the sum W Gn + W IGI + W IIG , and III to W In . The solid lines are evaluated with 
Az p = 0.0, the dashed lines with Az p = 0.1 and the dotted lines with Az p = 0.2. There is a non-negligible contribution to the sum 
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from the fraction p ca t of outlier galaxies at z = 0.5. 



2.2 Describing the redshift separation dependence 

The Az p dependence in the bispectrum as shown in Eq. [3] should primarily be determined by the weighting function W a ^ J , 
which has an explicit dependence on z p , z p , and z p . We find this to be particularly true of the ratio B(Az p )/B(Az p — 0). 
This observation forms the basis for understanding the Az p dependence in B GGG fSec. l2.2.1|) . B GGI (Sec. l2.2.4p . B Gn (Sec. 
I2.2.3|) . and B IU (Sec. 12.2.2"]) . We show in Fig. [3]the behavior of W aP ~' as a function of true redshift for various values of Az p . 

2.2.1 Redshift separation dependence in B 

We see from Fig.[3]that W , though strongly dependent on z p , has very little dependence on Az p , which explains the very 
small change in B as a function of Az p . In the 2-point case, Zhang describes the shear power spectrum with a Taylor 

expansion about Az p = up to second order, where the first derivative of the lensing spectrum with respect to Az p is shown 
to be zero. In the case of Method 3, it is not true that the first derivative is zero in the same way due to asymmetry. However, 
we can apply the same argument to our Method 1, which is symmetric, and express B as 



B 



GGG 



(Az p ) 



B GGG (Az p = 0) 

/GGG 



1 - f G GG(Az P Y, 

d 2 B GGG (Az p )/d(Az p f\ 
B GGG {Az p = 0) 



(13) 



We can apply the Jggg calculated in this way from Method 1 to compare with the results of Method 3. We find that for 
£ = 1000 and z p — 1, this results in Jggg = 0.67. Using this value, we find that the approximation in Eq. [13] is accurate to 
within 1% at Az p = 0.3. This is demonstrated in Fig. [4] where the right side of Eq. [13] is plotted as the dotted line labeled 
GGG, which nearly overlaps the left side, plotted as the solid line labeled GGG. In general, we find /ggg £ (0.5, TO) for 
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Figure 4. The ratio B a f 3 ~<(Az p )/B a P^(Az p = 0) evaluated for the labeled afif at I = 1000 and z p = 1.0. This demonstrates the 
accuracy of the relations in Eqs. 1131 & 1151 191 The lines labeled GGG which overlap are the left (solid) and right (dashed) sides of Eq. 
1131 The right side is evaluated as described in Sec. 12.2.11 The intrinsic alignment signal is evaluated using both the SB09 model and the 
toy modal with a = 1/2. For GGI+GIG+IGG and GII+IGI+IIG, the dashed lines use the SB09 model while the dotted lines use the 
toy model. The lines labeled III, Ilg, Igg, and ggg overlap due to the high accuracy of Eq. 1151 The relatively strong Az p dependences of 
the various bispectra compared to that of GGG can be understood by comparing the effects of different Az on the weighting functions 
in Fig. |31 which cause a strong decrease (or increase) in the amplitudes of W , W , and W IIJ . The relative differences in the 
lines associated with the GGI (GGg) and Gil (Ggg) bispectra are contributed to by both the different redshift dependencies of the two 
intrinsic alignment models, as well as the stronger dependence of the peak positions of W and W on redshift separation in Fig. 

[3] which decreases the accuracy of the approximation in Eqs. ll9l and ll8l 

£ £ (40,4000) at z p — 1. Alternatively we can fit Eq. [l3]to the calculated B and determine a best-fit fcGG- This results 

in a slightly lower value of /ggg with greatly improved accuracy, to within 0.1%. 



2.2.2 Redshift separation dependence in B 

Unlike W , it is clear from Fig. [3]that W 111 is very dependent on Az p . However, the peak position in z is determined by 
z p and is insensitive to Az p , as long as the photo-z is sufficiently accurate. The peak's amplitude decreases with Az p , which 
explains the quick decrease in magnitude of B 111 seen in Fig. [4] as Az p increases. This is a well-known result that is the basis 
for limiting measurements of the cosmic shear to only cross-correlations between thick photo-z bins in order to reduce the 
impact of the II, III, and Gil signals. 

This sharp peak in W 111 at z pca k, which is insensitive to Az p , allows us to make an approximation in the expression for 



B 1 



rjHI/n P P P\ Bin 

B {£; z 1 ,Z2,z 3 ) = 



(fc;X(2pcak)) 



W IU (x';Xi,X2,X3)dx, 



(14) 



X 4 (2pcak) 

where the dependence on Az p is contained within the final integral. In the limit of perfect photo-z information, this Az p 
dependence is exact. It is now clear that since W 111 = W 999 , we can relate B 111 and B 999 as 

B UI (£;z[,Z2 ,z£ ) = Aui(£;zi ,z£ , zf )B 999 (^; zf , z[ ,z£ ), (15) 

where Am = Biii(k;x( z paLk))/B g gg(k;x( z pea,k)). We propose to work with this relationship because it is in practice more 
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accurate than Eq. 1141 and because the same lensing survey will measure B 333 in addition to the shear bispectrum. This tells 
us directly the Az p dependence in B 111 without any knowledge of the photo-z PDF or intrinsic alignment being necessary. 
In the same way, we can construct the following relationships as well: 

B Il3 (£- z p , zf, 4") = A Ilg (£- z[, z p , z p )B 333 (e- z[ , z? , z p ), (16) 

B l33 (l; z P , 4, z P ) = A Igg (£; z p , z P , z p )B 333 (£; z[ , z P , z p ). (17) 

Figure [4] compares the Az p dependence of B UI , B IIg , B Igg , and B 333 and shows that Eqs. 1151171 are accurate to within 1% 
at Az p = 0.2. As with Jggg, we will treat Am, Au g , and Ai gg as free parameters in order to avoid modeling uncertainty 
in the self-calibration discussed in Sec. [3] 



2.2.3 Redshift separation dependence in B 

We find that the sum W + W + W 11 behaves very similarly to W 111 , as shown in Fig. [3] except for a slightly greater 
dependence in peak position due to Az p and much less dependence on Az p in peak amplitude. We thus propose a similar 
relationship to that described in Eq. [15] for B 111 , which is further motivated by the Az p dependence of B GI1 + B IGI + B IIG 
in Fig. [4] when compared to B " + B 9 3 + B" . Leaving out explicit dependence on £ and z p , z P , z P , this relationship is 
approximately 

B GII + B IIG + B IIG ^ Agii \ B Gg a + B gG a + ^gGl _ (lg) 

Due to the peak being less sharp as compared to that for B 111 and there being a stronger dependence on Az p in the peak 
position, this relationship is not as exact as Eqs. 1151171 We still find a large degree of accuracy, however, with Eq. 1181 being 
accurate to within 5% at Az p = 0.2. 



2.2.4 Redshift separation dependence in B 

Unlike W IU and W GU + W IGI + W IIG , the peak of the sum W GGI + W GIG + W IGG does clearly depend on Az p as shown 
in Fig. [3] though with a greater dependence in peak amplitude on Az p compared to the sum W + W 1 + W 11 and 
W 111 , instead increasing in amplitude with decreasing Az p . There is also a non-negligible contribution from the fraction 
Peat of outlier galaxies in Eq. [7] around 2 = 0.5. However, from Fig. 21 we still find a similar dependence on Az p between 
b ggi + b gig + b igg and B GGg + B GgG + B gGG ^ SQ we propose tne relationship 

B GGI + B GIG + B ,GG ^ Aqgi \^GGg + B GgG + B gGGJ (lg) 



While this relationship is much less accurate than Eos. 11511181 due to the assumption that z pe ak is constant being less valid for 
GGI, it still is accurate to withi n 5% at Az p = 0.1 and 20% at Az p — 0.2. Thus if we wish to work at the accuracy possible 



for the 2-point case as shown by I Zhang I (|2010bl ). where at Az = 0.2 the largest inaccuracy is about 10% in the GI term, we 
must limit ourselves to Az p ~ 0.13. Alternatively, we can use the limit proposed by Zhang of Az p w 0.2, but with a loss of 
accuracy in Eq. [19] of a factor of 2 when compared to the 2-point GI-Gg relationship. 

We can see from Fig. [4] that the accuracy of Eq. [19] could be improved by making assumptions about the redshift 
dependency of the intrinsic alignment or by measuring it independently. We could also include in the coefficient Aqgi an 
extra dependency on redshift, instead of the simple scaling relationship we have assumed, which could account for redshift 
dependent sources of error like the changing position of z pea k for GGI. This would introduce more free parameters into the 
self-calibration but would also allow the use of larger Az p . The choice of whether to include either intrinsic alignment model 
assumptions, additional measurements of redshift dependencies in the intrinsic alignment, additional free parameters in the 
self-calibration, or some combination of these or other assumptions will impact both the choice of upper Az p and the ultimate 
performance of the self-calibration, but this choice must be made with the requirements of a particular lensing survey in mind. 
For this reason, we will limit our discussion to the basic framework developed in this section when addressing the performance 
of the self-calibration in Sec. 13.11 



3 INTRINSIC ALIGNMENT SELF-CALIBRATION 

We can now consider a means to self-calibrate the intrinsic alignment signal, by using the four measurable bi spectra between 



galaxy ellipticity and galaxy density in a weak lensing survey, with negligible magnification bias, as discussed in lTroxel fc Ishak 
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Figure 5. The ratio B ( - 1 1(Az p )/B ( - 1 '> (Az p = 0), where £W = B GGG + B GGI + B GU + B In , and the intrinsic alignment signal is 
modeled with both the SB09 model (dashed) and toy model with a = 1/2 (dotted) for I = 1000 and z p = 1.0. From top to bottom, 
levels of total intrinsic alignment contamination |BW — B GGG \/B^ of 1, 5, 10, and 20% are shown. We are interested in the total 
absolute impact of the intrinsic alignment as a fraction of i?' 1 ), since this (and not the individual intrinsic alignment components) is 
what impacts cosmological information. For comparison we show B GGG (solid) with the expected minimum measurement uncertainty 
for an LSST-like survey. At this level of uncertainty, we find it is possible to identify the presence of even a small intrinsic alignment 
contamination of a few percent at Az p = 0.2. 



( 2012al ): 
B« = B 

B (2) = B GGg + g perm + B GIg + g p@mi + ^Ilg 



ggg + b ggi + 2 perm + B au + 2 perm + B ui ^ 



{;.',) 



(4) 



B G " +2 perm. + 3B 7 ", 
B aaa 



(20) 
(21) 
(22) 
(23) 



Equation [20] is the measured galaxy ellipticity-ellipticity-ellipticity bispectrum, which measures both cosmic shear and cor- 
related intrinsic alignment. In the case of no intrinsic alignment contamination, B^ 1 ' is simply the shear bispectrum, B , 
and should be eflectively independent of Az p within measurement error. However, we expect S' ' to have a measurable 
dependence on Az p at minimum survey error with even just a few percent or more contamination by intrinsic alignment due 
to the very strong Az p dependence of the intrinsic alignment bispectra. We explore this in Fig. [SI where we plot both B GGG 
and B^ for 1, 5, 10, and 20% levels of intrinsic alignment contamination for both the SB09 and toy model. We are interested 
in the total absolute impact of the intrinsic alignment as a fraction of B' 1 ', since this (and not the individual intrinsic align- 



ment components) is what impacts cosmological information, and so we measure the contamination as \B 



(i) 



B c 



; |/B (1) . 



Error bars representing the expected minimum measurement uncertainty in an LSS T-like survey a re shown on B . The 
measurement uncertainty in B is extrapolated from the e xpected error found bv lZhang I (2010bJ) for the power spectrum 

of LSST, modifying the derivation of iTroxel fc Ishak I (|2012al ) for the lensing bispectrum in a single redshift bin. 

Given simultaneous measurements of Eqs. I20ll23l in small redshift bins of size ~ 0.01 at 6 or more Az , along with the 
relations in Eqs. ll3l and ll5lfT9l we can now simultaneously reconstruct B , B , B , and B 111 through the observables 
given in Eqs l2QI23l and with free parameters fGGG, Aggi, Agii, Am, Ai gg , and Au g . In practice, it is more useful to measure 
the total contamination B +B +B 111 , which can be determined to higher accuracy due to a partial degeneracy between 
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the three intrinsic alignment bispectra. As we show in Fig. [4] this is because all three cause B^ 1 ' to decrease with Az p . B 111 
and B +perm. clearly decrease with Az p , and though B + perm. increases with Az p , its magnitude is negative, and so 
it also causes B^ 1 ' to decrease with Az p . Though we have concentrated on a particular £ and z p , Az p has a weak dependence 
on scale and mean redshift, so that we expect the self-calibration to be generally applicable to a wide range of values. 

3.1 Performance of the self-calibration 

From Fig. [3] we see that to good approximation _B Q/37 samples the same cosmic volume for different Az p , which means that 
different Az p should share the same cosmic variance as found for the power spectrum ( Zhang Il2010bl ). Relative differences 
due to cosmic variance are thus < V^i^Al 3 /^)-^ 2 = 1.4 x 10- 5 %(10 3 /^) 3/2 (i07 A ^) 3/2 / s ky /2 > and so are negligible 
compared to the > 5% change in B"' over Az p = 0.2 for intrinsic alignment contaminations of > 1%. However, the shot 
noise due to random galaxy shapes is large compared to the co smic variance at l arge £, for small redshift bins ~ 0.01 with 
a typical number of galaxies 2.5 x 10 7 at z p = 1.0 for LSST ( Zhan et al. 1 12003 ) . It must be controlled in order to reach 



the necessary precision in order to identify the Az dependence at low levels of intrinsic align ment contamina tion. This is 
complicated by the fact that unlike cosmic variance, shot noise is uncorrelated at different Az p .[ZhangJ (2010bJ) discusses a 
means to reduce the shot noise by averaging over larger £ bins at varying z p , since the Az p dependence is weak across £ and 
z p for the power spectrum. We find this to be accurate for the bispectrum as well, and expect this method to be applicable 
in the 3-point self-calibration in order to reduce the shot noise to manageable levels in the small photo-z bins necessary for 
the self-calibration. We also anticipate the self-calibration to be applicable in the presence of other errors which have different 
Az dependencies. 

In order for the self-calibration to reconstruct the intrinsic alignment signal in this way, with some associated error which 
is due primarily to the inaccuracy in Eqs. [19] & 1181 we must be able to measure 6 or more Az p in ~ 0.01 photo-z bins 
which are distinct from the measurement error. If we consider the expected measurement error for LSST shown in Fig. [5] we 
would be able to achieve the necessary measurements to self-calibrate the intrinsic alignment signal for intrinsic alignment 
contaminations of ~ 10% or more for maximum Az p = 0.2. If we instead accept a greater inaccuracy in Eqs. [T5]&[TH]at higher 
Az p , it would then be possible to measure smaller intrinsic alignment contaminations. The precise choice between maximum 
Az p and error in the final intrinsic alignment measurement through the self-calibration will then be entirely dependent upon 
the specific capabilities of the survey as well as the goals of the measurement. For example, if the resulting intrinsic alignment 
bispectrum is intended to constrain models of structure formation or intrinsic alignment models, then a lower maximum Az p 
might be imposed in order to achieve better accuracy in the intrinsic alignment measurement. 



3.2 Other sources of uncertainty 

The results of the previous section are dependent upon the accuracy of our assumptions in the quantitative calculations 
regarding the performance of the self-calibration. For the bisp ectrum, we have used an approximate fitting formula derived 
from perturbation theory bv lScoccimarro fc Couchman I (2001) in our performance estimations, and we have chosen a specific 
set of intrinsic alignment models. We also use a deterministic approach to modeling the galaxy bias, which is not perfectly 
accurate in real galaxy distributions. Failures in either these assumptions or the accuracy of the models for the intrinsic 
alignment bispectra would lead to additional uncertainty in the expected performance of the s elf-calibration, which may 
impact its applicability depending on the degree to which the assumptions or models fail. However. iBaldauf et al. I ( 20101 ) has 
shown that it is possible to suppress the galaxy stochasticity to the 1% level in some cases, which would be safely negligible 
compared to other sources of error we have discussed above. We have also chosen two very different models of the intrinsic 
alignment, in order to minimize bias in the resulting performance evaluation. 



4 CONCLUSION 

The strong Az p dependency of the intrinsic alignment signal for large bin size ^ 0.2 has been used to motivate the preference 
of cross-spectra and bispectra between redshift bins in order to reduce the intrinsic alignment contamination. This has 
previously been used as a means to neglect the III and Gil bispectra and II spectrum in techniqu es developed to remove 



intrinsic alignment contamination Irom the Icnsmg signal using mlormation between redstntt bins (Joachimi & achneidcr 
1, 120091 . l2010l : IShi. Joachimi fc Schneider II2OI0I : Izhang lEoiOaL iTroxel fc Ishak Il2012al ). In this work, we instead use this 



Az dependency within a single redshift bin of size ^ 0.2 to self-calibrate the measured galaxy ellipticity-ellipticity-ellipticity 
bispectrum, reconstructing simultaneously not only the cosmic shear bispectrum (GGG), but also the intrinsic alignment 
GGI, GII, and III bispectra. 

We first explore several means of defining and measuring the Az p dependence in the lensing and intrinsic alignment 
bispectra B , B , B , and B 111 , and show in Fig. [4] the resulting Az p dependency of the best method as well as 
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the accuracy of the scaling relations in Eqs. 11541191 These relate the intrinsic alignment to the directly measurable galaxy 
density through a set of simple scaling parameters. This is true for both the toy model and the SB09 model of intrinsic 
alignment clustering, which have very different dependencies on both scale and redshift. We further show that unlike the 
intrinsic alignment bispectra, B is effectively independent of Az p , which means a measured Az p dependence in the 

measured ellipticity bispectrum is a clear indication of the presence of intrinsic alignment. 

Using the measured cross-correlation galaxy ellipticity and galaxy density bispectra in a weak lensing survey, as well 
as the scaling relations in Eqs. I15H191 we propose a simple self-calibration method to simultaneously reconstruct the GGG, 
GGI, Gil, and III bispectra. The proposed self-calibration method relies only on the information already gathered by a weak 
lensing survey and makes no assumptions on the modeling of intrinsic alignment or photo-z PDF, instead depending on the 
relationship between the intrinsic alignment and galaxy density signals. We explore the feasibility and proposed structure 
of the self-calibration, including limitations on maximum Az p due to the inaccuracy of Eqs. [19] & 1181 and discuss some 
means of addressing it. We find that Az p ^ 0.2 is sufficient for accuracy in Eq. [19] better than 20%. The precision of survey 
measurements of the lensing bispectrum also limits the applicability of the self-calibration, but we find that for the expected 
measurement error in LSST, we can reconstruct the intrinsic alignment signal for contaminations of 10% at Az p — 0.2, or 
for even smaller intrinsic alignment contaminations, but with less accuracy in the reconstructed intrinsic alignment signal. 
Thus this would allow the self-calibration technique to significantly reduce the contamination of the intrinsic alignment to 
the weak lensing spectrum and bispectrum. However, further work is still necessary to precisely evaluate the quantitative 
performance of the self-calibration method in a realistic survey and its anticipated reduction in the effects of the intrinsic 
alignment contamination on cosmological study. 

The proposed self-calibration is complimentary to existing proposals for estimating the intrinsic alignment contamination 
to the bispectrum, which instead depend on information between redshift bins. It can be combined with these other methods 
which make use of redshift tomography to better constrain the intrinsic alignment contamination without over-using the 
information contained within the survey. By allowing the full reconstruction of the GII and III bispectra in addition to 
the GGG and GGI bispectra, while not relying on assumptions of an intrinsic alignment model, it also presents a means 
of indirectly measuring the intrinsic alignment signal, which will be applicable to constraining proposed intrinsic alignment 
models and our understanding of structure formation. 
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